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ABSTRACT 

Specification  of  parameters  defining  the  mean,  as  well  as  the 
turbulent  properties,  of  the  over  water  surface  layer  is  examined.    A 
baroclinic  boundary  layer  model  is  considered  using  Cardone's  modifi- 
cations for  the  marine  environment. 

Data  fields  of  pressure  and  temperature  were  selected  as  input, 
and  original  ship  observations  were  utilized  to  verify  the  predicted 
surface  parameters.    The  effects  of  stability  changes,  isobaric  curva- 
ture, thermal  wind,  and  the  value  of  the  von  Karman  constant  in  the 
model  were  examined  with  results  shown. 

Iteration  criteria  used  in  the  model  were  evaluated  with  respect 
to  operational  application.    Winds  calculated  by  the  model  at  the  height 
of  19.5  meters  were  incorporated  into  the  current  FNWC  wave  height 
program,  and  the  results  compared  with  observed  wave  data.    Turbulent 
parameters  describing  the  properties  of  the  index  of  refraction   (Cr)  were 
also  computed  and  evaluated  with  respect  to  laser  propagation  criteria. 
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I.      INTRODUCTION 

Parameterization  of  the  boundary  layer  over  the  sea  surface  has 
been  receiving  increased  attention  in  recent  years.    Fortunately  much 
of  the  theory  developed  for  low  level  turbulence  over  a  land  surface  is 
also  applicable  over  the  oceans.     But  the  air-water  interface  does  pre- 
sent some  unique  problems.    For  instance  a  fluid  surface,  rather  than  a 
fixed  one,  results  in  mobile  surface  roughness  elements.     Because  of 
the  presence  of  a  fluid  surface,  any  parameterization  requires  a  knowl- 
edge of  the  processes  which  transfer  heat  and  momentum  energy  across 
the  boundary  and  ultimately  influence  such  variables  as  wave  generation 
and  index  of  refraction. 

The  purpose  of  this  study  was  to  examine  and  evaluate  several 
features  of  a  recently  formulated  marine  boundary  layer  model  which  was 
designed  for  operational  use.    The  model  is  a  two-layer  baroclinic  model 
for  the  marine  boundary  layer  and  the  formulation  examined  was  that 
developed  by  Cardone  (1969).     Cardone  based  his  formulation  on  an 
earlier  two-layer  model  due  to  Blackadar  (1965a).     Features  of  interest 
in  the  model  are  related  to  the  physics  of  the  planetary  boundary  layer 
and  include  the  effect  of  atmospheric  stability,  baroclinicity,  coriolis 
acceleration  and  changes  in  the  surface  roughness  parameter  on  the 
pressure  driven  surface  wind. 
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The  scope  of  the  study  was  to  examine  the  model  in  the  context  of 
its  operational  usage.    As  such,  it  was  a  case  study  in  which  the  data 
represented  pressure  and  temperature  observations  from  sequential 
12-hour  periods  when  a  large  cyclonic  pressure  system  dominated  the 
Central  North  Pacific  Ocean.    This  particular  synoptic  situation  provided 
the  necessary  relations  between  the  scalar  fields  to  examine  features 
associated  with  baroclinicity,  stability  and  wind  speed  categories. 

Because  these  data  represent  scalar  fields  which  had  to  be  de- 
fined by  some  analysis  scheme  before  being  used  as  input  to  the  two- 
layer  model,  the  differences  which  arise  due  to  procedures  in  the  external 
field  specifications  are  considered.    The  effect  of  isobaric  curvature  on 
the  surface  wind  fields  is  examined  and  this  effect  is  compared  to  the 
changes  associated  with  the  physics  in  the  model.    The  differences 
which  may  arise  solely  because  the  fields  are  defined  by  an  objective 
versus  a  subjective  analysis  are  considered. 

Finally,  expressions  describing  the  boundary  layer  are,  by  their 
very  nature,  empirical.    As  such,  they  contain  empirical  constants   and 
the  von  Karman  constant  (k)  is  one  whose  value  has  recently  been  ques- 
tioned.   The  importance  of  a  suggested  change  from  k=0.40  to  0.35  was 
considered. 
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n.      ASPECTS  OF  A  PLANETARY  BOUNDARY  LAYER 

A.         SURFACE  BOUNDARY  LAYER 

1.         Basic  Structure  of  the  Boundary  Layer 

Recent  advances  in  the  study  of  atmospheric  turbulence 
have  led  to  a  greater  confidence  in  describing  the  distribution  of  wind 
and  temperature  in  the  atmospheric  boundary  layer  over  a  solid  surface. 
This  progress  has  arisen  largely  through  the  application  of  the  Monin- 
Obukhov  similarity  theory  which  has  been  supported  by  recent  observa- 
tional studies.     Cardone's  contribution  was  primarily  to  extend 
Blackadar's  two-layer  neutral  baroclinic  boundary  layer  model  over  fixed 
terrain,  and  then  to  apply  this  model  to  the  boundary  layer  over  a  fluid 
surface.    The  major  difference  between  these  boundaries  is  that  the  sea 
surface  changes  its  character  as  the  wind  blows  over  it. 

In  general,  Blackadar's  model  treats  the  boundary  layer  as 
being  formed  by  both  turbulent  friction  and  the  coriolis  acceleration.  All 
processes  are  considered  to  be  three  dimensional,  and  the  layer  is  char- 
acterized by  a  density  stratification.    It  is  conveniently  separated  into 
the  three  regions  depicted  in  Figure  1 .       In  this  delineation  the  important 
consideration  is  the  specification  of  Km  which  is  the  turbulent  transfer 
coefficient  for  momentum. 
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Fig.   1.      Diagram  of  the  planetary  boundary  layer 

a.         Viscous  Sublayer  (less  than  1  cm) 

In  this  near  surface  layer  the  stress  is  supported  by 
viscosity  (aerodynamically  smooth  flow),  and  the  mean  velocity 
increases  linearly  according  to 


* 


u    = 
z 


(1) 


where  v    is  the  kinematic  viscosity  and  U^  is  the  friction  velocity. 

b.         Surface  Layer  (up  to  50  meters) 

This  is  the  constant  stress  layer  in  which  the  coriolis 

parameter  can  be  neglected  and  K      increases  linearly  with  height.    For 

m 

neutral  stratification  the  logarithmic  law  is  valid 
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Uz   =     T   log     /zo  <2> 

where  k=0.4  is  the  von  Karman  constant. 

c.         Ekman  Layer  (50-600  meters) 

This  layer  is  characterized  by  the  stress  decreasing 

with  height  and  K     constant  throughout  the  layer, 
m 

In  formulating  the  boundary  layer  model,  Cardone  (1969) 
considered  the  validity  of  the  "log  law"  for  the  surface  layer  over  water, 
which  has  been  experimentally  documented  down  to  near  surface  limits 
(Phillips,   1966;  Paulson,   1967).     He  also  assumed  that  the  process  of 
wave  generation  did  not  seriously  affect  the  surface  layer  since  under 
neutral  and  steady  conditions  a  logarithmic  shear  zone  should  be  estab- 
lished.   Accepting  the  log  profile  as  being  valid,  Cardone  was  able  to 
apply  much  of  the  surface  layer  turbulence  theory  developed  over  land 
to  over-water  conditions .    However  Cardone  did  take  into  account  the 
action  of  the  waves  when  specifying  the  effective  roughness  parameter 
for  the  sea  surface. 

2 .         Similarity  Theory  for  the  Surface  Layer 

The  Monin-Obukhov  "  similarity  theory"  has  led  to  consider- 
able progress  in  the  specification  of  the  distribution  of  wind  and  tem- 
perature in  the  boundary  layer.    The  basis  of  this  theory  is  that  near 
the  ground  there  exist  velocity  (U^),  length  (L),  and  temperature  (Qj 
scales  which  are  essentially  invariant  with  height.    When  the  principle 
variables  such  as  temperature  (T) ,  wind  (U) ,  and  height  (z)  are  expressed 
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as  nondimensional  fractions  of  these  quantities,  a  series  of  nondimen- 
sional  equations  can  be  formed  which  are  of  general  validity  in  the  sur- 
face boundary  layer  (Lumley  and  Panofsky,   1964).    Recent  analyses  of 
the  oceanic  wind  data  (e.g.  Phelps  and  Pond,   1971)  have  indicated  that 
the  Monin-Obukhov  theory  is  valid  for  over-water  application.    These 
results  support  Cardone's  inclusion  of  the  Monin-Obukhov  stability 
criteria  in  the  boundary  layer  model. 

Similarity  theory  predicts  that  a  universal  relation  should 
exist  for  the  nondimensional  wind  shear 

U*        az    =    ^u  (3) 

and  the  nondimensional  temperature  gradient 

■a       I3"    =    *t  (4) 

8*       dz  x 

where  #t  and  0U  are  functions  of  the  dimensionless  ratio  z/L.    Here 

Q     is  the  potential  temperature,    0*  is  the  scaling  temperature  defined  by 


1        H_ 

kU*      Cp     Pa 


&*  =-T7T-     „■»_  (5) 


and  L  is  the  Lettau-Monin-Obukhov  stability  length  defined  by 

.     -Vl   Cp     Pa T 

kg  H 

The  heat  flux   (H)  is  not  easily  measured,  so  a  modified  stability  length 

was  defined  by  Panofsky  (1963) 

KH  U*T   du/dz 

L'  =  ahL  =  k^l  =  kg   ae/az  ■ 
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The  assumption  that  L'  is  independent  of  height  implies  that  a,    is  con- 

n 

stant,  which  in  turn  implies  similarity  between  wind  and  temperature 
profiles . 

If  L'  is  a  valid  scaling  length,  similarity  theory  predicts 


that 


fr   It  =  *«<**■>  (8) 

z      *0        Vz/L,) 


9*       dz  au  *  (9) 


Integration  of  these  formulae  (Panofsky,   1963)  yields 


U* 


Q 

9=9     +   — 
z        o  a 


In    z/zQ-  ^(z/L)  (io) 

—        In  z/zQ  -     ^(z/L1)      ,  (11) 


where  the  relation  between    H>    and  0     is 

/•Z/L'  i-0u(i) 

V(z/L')=/  r d*   ,      1=  z/L\  (12) 

•'o  * 

Paulson  (1970),  using  the  KEYPS  formula,  integrated  (12)  by  parts  to 
obtain  the  following  expression  for  *P 

V  =    l-J*u-31nj*u  +  21n   (ip)  +  2  tan"1^  -  -f.  +  ln  (^.(13) 

Within  the  scope  of  an  operational  application,  these 
expressions  for  the  surface  layer  would  be  used  to  determine  the  wind 
speed  given  U*  ,    or  to  determine  U*  given  the  wind  speed  at  a  particular 
level.    For  this  reason  the  partial  derivatives  in  (7)  are  replaced  by 
expressions   involving  height  differences   and  air-sea  temperature 
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differences.    This  is  achieved  by  substituting  the  following  expression 

d  9 

for       t —   ,  obtained  by  combining  (9)  and  (11) 
cz 

3z         z  [ln(zVz0)-*(Za/L')]   '  (14' 

where  zQ  is  the  height  at  which  8a  is  measured  (chosen  as  10  meters  by 

Cardone) . 

Substituting  (14)  and  (8)  into  (7)  yields 

z- 


uf  8   [in  (z§)-     ^(Za/L-)] 


l'  -  — — — * — : J  .  (is) 


k'g  (€a  -  9S) 

Equation  (10)  is  then  solved  for  U* 

k  um 

U*    = ffi ,  (16) 

ln(zm/z0)-     W(zTn/L>) 

where  z      is  the  height  at  which  the  wind  speed  is  measured  (chosen  as 
19.5  meters  by  Cardone).      Equations  (15)  and  (16)  can  then  be  solved 
for  U*  given  the  wind  speed  at  some  level  (U2)  and  the  air-sea  tempera- 
ture difference . 

3 .         Consideration  of  the  Drag  Coefficient  ( Cz)  and  Roughness 

Parameter  (zo) 

The  surface  roughness  parameter  and  the  influence  of  thermal 
stratification  have  received  considerable  attention  in  the  definition  of 
the  relation  between  the  nondimensional  drag  coefficient  (Cz)  and  the 
wind  speed  at  a  certain  height  (usually  10  meters).    The  need  for  such  a 
relation  arises  from  the  "bulk  aerodynamic"  formula  for  estimating 
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momentum  flux.    This  equation,  deduced  from  dimensional  analysis,  has 
the  form 

T  =     Pa    U*   =    Pa  Cz  Vl     <  <17) 


where    T   is  the  surface  stress.    From  this  expression  the  drag  coefficient 
is  defined  as 


u: 


c*=  1* 


(18) 


Rearranging  and  squaring  (10)  yields  the  following  expression  for  the 


drag  coefficient  at  a  level  z  in  terms  of  z     and  stability 


Cz    = 


In  Z/zQ-  vf/V) 


(19) 


which  was  the  form  considered  by  Cardone. 

Observations  have  indicated  that  over  the  sea,    zQ  is  depend- 
ent on  U*.      This  dependence,  however,  changes  as  the  wind  speed 
changes.    For  the  case  of  light  winds  (aerodynamically  smooth  flow), 
the  drag  coefficient  C     decreases  with  increasing  wind  speed  (Kraus, 


1966)  and  zn  has  the  form 


Zo  1.09   U*     ' 


(20) 


For  aerodynamically  rough  flow,    C     increases  with  increasing  wind 


speed  according  to 


TT2 
a  U* 


(21) 


where   (a)  is  a  constant  =  0.035  (Charnock,   1955).    A  summary  of  field 
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observations  indicates  that  the  change  between  rough  and  smooth  flow 

occurs  at  wind  speeds  of  5-6  meters/second. 

Cardone  proposed  the  following  simple  expression  for  z 

o 

covering  a  wide  range  of  wind  speeds  (both  rough  and  smooth  flow) 

Cl  2 

zQ   =    —      +     C2U*    +    C3  .  (22) 

The  constants  were  chosen  so  that  C      was  a  minimum  for  a  value  of 
6  meters/second  at  the  10  meter  height  under  neutral  conditions,  and 
above  this  speed  the  relation  approximated  Charnock's  results.    For  U* 
in  meters/second  and  z     in  meters,  the  resulting  expression  was 

zQ  =  6.84  x  10~5/TJ*  +  4.28  x  10_3U^  -  4.43x  10".4       (23) 

Other  representations  of  C     versus  U     have  been  published. 
A  few  of  these  appear  in  Figure  2  along  with  Cardone' s  approximation. 
It  is  clear  from  the  scatter  of  curves  in  Figure  2  that  the  specification 
of  the  drag  coefficient  was  a  primary  decision  in  Cardone' s  formulation. 
This  choice  of  C„   versus   U_   will  be  examined  during  the  discussion  on 
results  from  this  investigation. 

With  the  equations  for  stability  length  (15),  friction  velocity 
(16),  and  surface  roughness  parameter  (23),  the  surface  layer  wind  dis- 
tribution is  completely  described  within  the  framework  of  the  similarity 
theory.    These  three  equations  can  be  solved  simultaneously  for  U* 
using  iteration  techniques  with  the  external  parameters  described  below 
as  input. 
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B.         PLANETARY  BOUNDARY  LAYER 
1.         Background 

In  developing  a  model  for  the  entire  boundary  layer,  con- 
sideration has  to  be  given  to  the  possible  external  parameters.    In 
Cardone's  formulation  these  parameters  were  the  geostrophic  wind  (G), 
the  air-sea  temperature  difference  (9    -  Q    ),  the  coriolis  parameter  (f), 
and  the  thermal  wind  (r) .    He  considered  previously  suggested  empirical 
relations  between  the  geostrophic  drag  coefficient   (U*/G)  and  the  sur- 
face Rossby  number  (Lettau,    1959) 

Ro   -    ff"         -  <24> 

1  ^o 

where  G  is  the  surface  geostrophic  wind  speed.    Empirical  relations 
have  also  been  suggested  between  the  cross  isobar  flow  angle  ip     and 

V 

Blackadar  (1965a)  developed  a  two- layer  boundary  layer 
model,  consisting  of  a  Prandtl  surface  layer  and  an  Ekman  layer,  which 
was  verified  on  the  basis  of  the  above  relations. 

Cardone's  formulation  was  closely  related  to  Blackadar' s 
work.  In  particular,  the  height  of  the  surface  layer  (h)  was  specified 
explicitly  in  terms  of  the  external  parameters 

h   =   ^  (25) 

-4 
where   BQ  was  a  nondimensional  constant  equal  to  3. Ox  10       as  deter- 
mined by  Blackadar. 
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By  applying  proper  boundary  conditions  and  requiring  that 
continuity  of  wind  and  wind  shear  exist  at  level  h,  the  following  solu- 
tions for  U*/G  and     tyQ   (along  with  the  solution  for  z    )  completely 
define  the  surface  boundary  layer  wind  distribution  under  neutral 

conditions . 

■,3/2 

(26) 


k  BQ  sin2     V0 


U^  k>/7     sin    (n/4-  ^q) 


(27) 


G  ln  Bo  Ro 

Implicit  in  such  an  approach  is  that  the  value  of  K      at  the  internal 
boundary  between  the  surface  layer  and  the  Ekman  layer  is  given  by 

kU*  G 
Km    =    — —  .  (28) 


2 .         Cardone's  Extension  of  Blackadar's  Two-Layer  Model 

The  success  of  Blackadar's  two-layer  model  in  predicting 
the  essential  characteristics  of  the  neutral  surface  boundary  layer,  along 
with  its  relative  ease  of  application,  suggested  its  extension  to  the  non- 
neutral  marine  boundary  layer.     Cardone  extended  the  diabatic  model  to 
over- water  flow  by  regarding  zQ  as  a  function  of  U*  . 

The  eddy  viscosity  in  the  constant  stress  layer  for  the  non- 
neutral  case  was  expressed  by  Cardone  in  terms  of  the  nondimensional 

wind  shear  as 

kU*z 

^   = •  (29) 

0U(ZA') 
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Since  Km  is  constant  with  height  above  h,  its  value  in  the  Ekman  layer 

is  then 

kU.B^G 

(30) 


K      =   l_!f 

»         f0u(z/L') 


The  wind  speed  at  h  was  given  by 


h  k 


lnBQRo-    m  (h/L1) 


(31) 


Internal  boundary  conditions  applied  at  level  h  appear  in  Figure  3  and 

r  ,,  ,  ,  dw      ,  ,  3n 

follow  from  the  fact  that  always  makes  an  angle  of     — - —    to 


dz 

the  ageostrophic  wind    (W)  and  that 


constant  stress  layer. 


is  parallel  to  U  in  the 


Fig.  3.     V/ind  vector  diagram  applicable  at  level  h 


By  requiring  continuity  of  wind  shear  at  h,    as  well  as  directional 

continuity,  Cardone  was  able  to  reduce  the  problem  (analogous  to  (26) 

and  (27))    to 

1/2 


-^=     [2  k  BQ  sin2    V0£u(hA')] 

U*  _       kyJT  sin    (  n/4  -M*0) 
G  In  B0R0  -   i|i  (h/L1) 


(32) 
(33) 
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These  equations  are  solved  simultaneously  for   —   and    ip        from  the 

G  o 

external  parameters  G,    f,    and   (Qa~Qs)«      In  addition,  a  dimensionless 
stability  parameter  was  defined  as 

B0G 

L*  =  TIT  (34) 

which  represented  the  ratio  of  the  surface  layer  h  to  the  stability 
length  L'  (Blackadar  and  Ching,   1965). 

The  relationships  determined  by  Cardone  for  U*/G  and  ip 
as  functions  of  RQ  and   L*  appear  in  Figures   4  and   5.      For  all  surface 
Rossby  numbers ,    U*/G  increases  and      ID      decreases  with  decreasing 
stability.    In  general,  stable  conditions  have  a  more  pronounced  effect 
on  the  surface  layer  wind  characteristics  than  unstable  conditions.    In 
comparing  the  results  from  the  model  with  observations,  Cardone  noted 
that  U*/G  increased  by  roughly  30%  in  unstable  conditions,  with  no 
marked  tendency  to  decrease  at  great  instability,  and  decreased  roughly 
70%  in  stable  conditions.    These  features  in  the  observed  results  will 
be  discussed  later. 

The  effects  of  baroclinicity  were  also  included  by  assuming 
that  the  geostrophic  wind  was  a  linear  function  of  height  (Blackadar, 
1965b).     Blackadar  (1965a)  introduced  a  parameter,  p,  which  was  defined 
at  level  h  by 


Uh  BoG 

p  = 


(du/dz)h      fjefu(h/L') 


In  B0Rq  -iP(hA')|  .       (35) 


'] 
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Fig.  4.  Variation  of  the  geostrophlc  drag  coefficient 
with  surface  Rossby  number  and  dimensionless 
stability  length  (1^)  (Cardone,  1 969) 
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The  requirement  of  continuity  of  wind  shear  for  the  baroclinic  case  took 
the  form 

(36) 


dW 


lf)+ 


1- 


Figure  6  illustrates  the  geometrical  implication  of  (36).    From  the  Ekman 
layer  solution  of  the  ageostrophic  wind,  the  following  equation  can  be 
obtained 


dw 


V2*  p'    I  wh  I      , 


(37) 


where  p1    is  a  dimensionless  parameter  defined  by 

0     (VL1)  r 


P'   = 


B    G 
o 


0u(h/L>) 


{ 


u 


In  B0R     -   W  ( 


2kU+B    G  °    o 


h/L' )  1 


(38) 


aw 


Fig.  6.      Geometrical  relationships  resulting  from  the 
requirement  of  the  continuity  of  wind  shear  at. 
the  internal  boundary  level  h  (Blackadar,  1965b) 


Next  the  thermal  wind  vector  was  replaced  by  its  dimensionless 
magnitude 


r    = 


f 


ac 


(39) 
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where  the  direction  of 


dG 


relative  to  the  surface  geostrophic  wind 


was  implied  by  angle  n    in  Figure  6.      Still  another  dimensionless 
quantity  was  defined 


r'    =    rB      r  In  B    R     -   IV  (h/L')l 
o    L  o    o  J 


(40) 


and  (35)  was  then  rewritten  as 


dG 


=     r' G 


(41) 


By  applying  the  law  of  cosines  to  each  of  the  two  smaller  triangles  in 
Figure  6  that  have   M  for  one  side,  the  following  relations  were  determined 

2 


M      =        Wi 


2  -2    2 

s      =    G    q 


s2     =      (1  +  2p'  +  2p'2) 


=      (1  +  r'2  -  2r'  cos  H  )• 


(42) 
(43) 
(44) 


H    is  the  angle  between  the  geostrophic  wind  and  the  thermal  wind  (see 
Figure  6).     Then  the  law  of  sines  was  applied  to  yield 

sin  d  P'/s  (45) 

sin  7        =       —  sin  f  (46) 


—  sin  n 

q 


The  following  relations  apply  to  the  large  triangle  (ABC)  in  Figure  6. 


Uu  =     |Wu       sin    (  a  +  7  )/sin  Ui 


U2=    G2  + 


Wh       -2G     Wh|  cos  (  a   +  7  ) 


W, 


2    =    U?+G2-2UuGcoslii 
n  ho 


(47) 
(48) 
(49) 


From  equations  (31),   (45),   (46),   (47),   (48),  and    (49),    the   following 
equations  describing  the  surface  boundary  layer  characteristics  were 


obtained 
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sin  (  a    +   7    ) 

Tan   m     =    : : (50) 

0  (s/q)-  cos  (  a  +  7  ) 

U      =        ^  G  q  sin  (  a     +   7    ) 
*  s    sin  Ui     [In  BR-   Ul(h/L')]    '  K     ' 

By  including  (23)  and  (15),  Cardone  obtained  a  system  of  equations  that 
could  be  solved  for  U*  and    HJ  _  from  the  input  parameters  G,  f, 


(e  -©J, 


o 
IdG 


,    and    r|   .       These  equations  were  calculated  from  a 


a       s ' '        fez 
knowledge  of  sea  level  pressure,    surface  air  temperature,  sea  surface 

temperature,  and  latitude. 

3 .         Some  Characteristics  of  the  Marine  Surface  Boundary  Layer 
A  primary  effort  of  Cardone' s  work  was  to  determine  how 
well  the  planetary  boundary  layer  model,  using  only  large-scale  synoptic 
parameters  as  input,  specified  the  surface  boundary  layer  wind  distribu- 
tion over  the  sea  surface.     Cardone  considered  several  studies    relating 
the  geostrophic  wind  ratio  Uz/G  and  the  air-sea  temperature  difference. 
The  results  of  an  extensive  study  by  Carstensen  (1967)  at  Fleet  Numerical 
Weather  Central  (FNWC)  appear  in  Figure  7  along  with  earlier  results 
due  to  Bleeker.    The  ratio  of  surface  wind    (V)  to  geostrophic  wind   (G) 
is  tabulated  for  each  degree  of  the  reported  air-sea  temperature  differ- 
ence, and  the  median  value  for  each  interval  is  plotted  as  an  x  in 
Figure  7.    The  validity  of  the  results  was  limited  to  the  air-sea  tempera- 
ture difference  range   -4   to  +1°C.     In  order  to  include  observed  winds 
in  the  surface  pressure  analysis,  FNWC  uses  a  constant  value  for  V/G 
of  0.78  and  a  turning  angle  of  15°  toward  lower  pressure.    This  selection 

could  produce  significant  error. 
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Cardone's  model  can  be  used  to  compute  a  wind  speed  at 
19.5  meters.     Figure  8  illustrates  the  variation  of  the  ratio  of  the  19.5 
meter  wind  and  surface  geostrophic  wind  with  the  air-sea  temperature 
difference.    The  ratio  varied  less  in  unstable  conditions  than  in  stable 
conditions . 

The  effects  of  baroclinicity  also  appear  in  Figure  8  and 

these  effects  increase  as  stability  decreases,  with  the  ratio,  U         /G , 

19  . 5 

increasing  with  cold  advection  and  decreasing  with  warm  advection. 
Mendenhall  (1967)  also  studied  the  influence  of  baroclinicity  and  sta- 
bility on  the  veering  of  the  wind  in  the  planetary  layer.    He  concluded 
that  over  mid-latitude  oceans  the  diurnal  variation  of  frictional  veering 
of  the  wind  with  height  was  negligible,  but  that  lapse  rate  and  baro- 
clinicity were  important.    Results  obtained  by  Cardone,  using  the  plan- 
etary model,  showed  qualitative  agreement  with  Mendenhall1  s  conclu- 
sions.   Cardone's  results  also  suggested  that  baroclinicity  was  more 
significant  than  stability  in  determining  the  veering  of  the  wind  in  the 
planetary  boundary  layer.    For  both  cases  in  situations  of  strong,  cold 
advection  the  actual  wind  may  back  with  height,  especially  in  conjunc- 
tion with  unstable  stratification. 

Additional  results  on  the  influence  of  baroclinicity  came 
from  a  study  by  Clarke  (1970),  who  found  a  tendency  for  the  boundary 
layer  wind  to  back  instead  of  veer  with  height  in  unstable  conditions . 

He  stated  that  this  could  have  been  due  to  baroclinicity  or  to  accelera- 

3  T* 

tion.    He  noted  that  it  could  have  occurred  if     ~^r        was  systematically 

o  z 
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u* 

negative  in  convective  conditions  up  to  a  level  of  about  0.2   ~7 — 
Instead  of  being  positive  as  in  the  normal  case.    This  latter  interpreta- 
tion was  based  on  the  observation  that  the  wind  appeared  to  be  sub- 
geostrophic  in  the  convective  layer  and  reached  geostrophic  values  only 
after  a  rapid  increase  through  the  overlying  inversion. 
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III.      COMPUTATIONS  AND  NUMERICAL  PROCEDURES 

A.         EXTERNAL  WIND  FIELD  SPECIFICATIONS 

All  computations  were  made  with  data  defined  on  a  subset  of  the 
FNWC  63X63  grid    which  covered  the  North  Pacific  Ocean  (Figure  9). 
Data  fields  included  sea  level  pressure,  sea  surface  temperature,  and 
surface  air  temperature  at  each  of  315  grid  points  spaced  381  km  apart 
at  60°N.    The  necessary  map  factor  was  included  for  a  1:30  million 
polar  stereographic  projection. 

The  coriolis  term  appears  as  an  external  parameter  and  therefore 
poses  a  limitation  on  the  lowest  geographical  latitude  where  the  geo- 
strophic  relation  would  be  valid.    It  has  already  been  established  that 
the  dynamic  boundary  layer  was  formed  by  the  surface  stress  and  the 
coriolis  acceleration.    Pushistov  (1970)  concluded  that  manifestations 
of  the  nonlinearity  in  the  dynamics  of  the  boundary  layer  and  the  appear- 
ance of  a  perceptible  ageostrophic  component  of  velocity  in  the  free 
atmosphere  were  to  be  expected  near  latitudes  of  7-10°.      Sheppard 
(  1970)  suggested  20°N  as  a  limiting  latitude  for  application  to  boundary 
layer  models.    A  limiting  value  of  15°  latitude  was  used  in  this  program 
as  the  limit  on  geostrophic  flow.     For  latitudes  below  15°,  the  coriolis 
parameter  would  be  kept  constant. 

1.         Geostrophic  Wind  Computations 

The  geostrophic  wind  components  were  computed  from  the 
surface  pressure  field   (P)   as  follows 
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u      = 

g 

paf 

ay 

1 

dp 

V        — 

g 

P* 

dx 

dP  /     x 

(52) 


(53) 


using  a  centered  explicit  finite  difference  technique.    The  density  was 
determined  from  the   equation  of  state 

Pa   =    R-^  <54> 

where  R ,   is  the  gas  constant  for  dry  air,  and   0    the  surface  air  tempera- 
ture in  degrees  Kelvin.    The  magnitude  of  the    geostrophic   wind   was 
computed  at  each  grid  point  from 


V 


f  u2    +    v2  .  (55) 

g         g 


2.         Correction  for  Curvature 

A  gradient  wind  was  calculated  to  correct  for  the  effect  of 
curvature  on  the  surface  geostrophic  wind.  Two  basic  principles  were 
required  for  this  correction.  The  first  was  that  following  an  air  parcel 
along  an  isobar,  the  total  derivative  was  zero  or 

dP(x,y)    =    0.  (56) 

The  second  requirement  was  the  mathematical  statement  of  curvature 

Kc(x'y)  =     [l  +  (y/]V2     °i  (57) 

where  y'  and  y"  were  the  first  and  second  derivatives  with  respect  to  x, 
and  R  was  the  radius  of  curvature.    By  expanding  the  total  derivative  of 
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P,  obtaining  y'  and  y" ,  and  rearranging,  the  following  result  was  ob- 
tained for  the  radius  of  curvature 


3, 

'2 

sp    [„  a2  p         d2  p 

d  x            dx  dy            Sy^ 

c)P/c-x~ 
aP/c3y 

dP    d2P 
c3y      dx2 

R  =  r — r? — — ~5 r— n — : — ^r~  •    (58) 


Next  the  gradient  wind  was  computed  from 

V2 
VG    =    G  -    _G_     .  (59) 

fR 

Solving  for  V      and  rearranging  led  to  a  workable  form  of  the  equation 

v,~    =    i  /o  j. — i  =■     .  (60) 

G  V2  +  ^  1/4  +  G/Rf 

For  this  case  the  magnitude  of  the  term  G/Rf  determined  how  much  curva- 
ture would  be  allowed.    The  maximum  values  were  chosen  to  be  -3/16 
for  anticyclonic  curvature  and  +3/4  for  cyclonic  curvature  which  yielded 
a  range  of  V\_,  =  1.33G   for  anticyclonic  curvature  and  V    =0.67G  for 
cyclonic  curvature.    This  agreed  with  generally  accepted  synoptic  limits. 

Subroutine  GRADWD  (Appendix  A)  makes  the  necessary  com- 
putations and  follows  a  similar  scheme  utilized  by  FNWC  (Kaitala,    1971 
unpublished).     Centered  finite  difference  schemes  were  used  to  calculate 
the  derivatives. 

3  .         Thermal  Wind  Correction 

The  nondimensional  thermal  wind  components  were  deter- 
mined from 

u      -    -       9  _  -—  (61) 

T  f  2  Q  ^y 

a 
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vt  =    fit"  Ti  '  (62) 


and  the  total  thermal  wind  was  given  by 


r    =  V  UT  +  VT      *  <63) 

For  this  model,    r  had  a  limiting  value  of  100.    The  geostrophic  and 
thermal  wind  angles  were  calculated  from  the  previously  computed  wind 
components .    This  procedure  was  repeated  for  each  grid  point  and  was 
used  as  input  data  to  the  boundary  layer  model. 

B.         PROCEDURES  IN  SOLVING  THE  BOUNDARY  LAYER  EQUATIONS 
1 .         Initial  Conditions 

a.         Geostrophic  Wind  Speed  Less  Than  Or  Equal  to  5 . 0 
meters/second: 

This  case  was  associated  with  the  aerodynamically 
smooth  flow  region  described  earlier.    For  this  case  no  corrections  were 
made  for  stability  and  a  constant  drag  coefficient  of  0.022  was  assumed. 
This  resulted  in  a  constant  inflow  angle  of  15°.    The  19.5  meter  wind 
was  determined  from 

U19.5    "    °'7G  (64) 

The  choice  of  the  drag  coefficient  was  not  arbitrary,  but  depended  on  the 
relationship  between  C     and  U     (see  Figure  2).    The  value  of  C   _   was 
a  minimum  at  6  meters/second  under  neutral  conditions.    The  value  of 
C1Q    _   was  extrapolated  from  this  curve. 

iy .  5 
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b.  Geostrophic  Wind  Speed  Greater  Than  5.0  meters/ 
second  and  Absolute  Value  of  Air-Sea  Temperature 
Difference  Less  Than  or  Equal  to  1°C: 

Initially  the  drag  coefficient  was  C    =  0.0245,  with 
the  thermal  wind   (r)  and  the  angle  between  the  geostrophic  and  thermal 
wind  (  H  )   set  equal  to  zero.    For  this  case  L*,  L'  and     ty  (h/L1)  were 
set  equal  to  zero  and  0u  (h/L')=  1.0.    A  solution  was  then  obtained  using 
the  modified  baroclinic  model  of  Blackadar  described  below. 

c.  Geostrophic  Wind  Speed  Greater  Than  5.0  meters/ 

second  and  Absolute  Value  of  Air-Sea  Temperature 

Difference  Greater  than  1°  C  : 

The  drag  coefficient  C    =0.0245  is  used  to  estimate 

z 

the  initial  U*  and  r  and    rj  are  set  equal  to  zero.      0-  -  Q     was  restricted 

a       s 

to  the  range  -15  to  +4.0 °C.    The  values  of  L+,  L' ,     ty  (h/L1),  and 

0     (h/L1)  were  computed  from  z    ,  0    -€     and  U     and  a  solution  was 
u  a       a       s 

obtained  using  the  modified  Blackadar  model  described  below. 

L^.   was  computed  from  (34)  and  L'  from  (15)  using  as 
a  first  approximation 

W(za/L')    =    10-  In  (zVz0)  .  (65) 

m(h/L')  was  computed  from  SUBROUTINE  PSI  (Appendix  A).     If  L^   was 
greater  than  zero  (stable  case),  then 

m(h/L')    =    0.7  L.  (66) 
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If  L^  was  less  than  or  equal  to  zero  (unstable  case),  a  value  of 
Richardson  number  and  0U  (h/V)  had  to  be  computed.    This  was  accom- 
plished in  SUBROUTINE  SHR  (Appendix  A).    Then    IJJ  (h/V)  was  calculated 
from  (13). 

0U  (h/V)  was  computed  using  SUBROUTINE  SHR 
(Panofsky,   1968;  Cardone,   1969).    If  L+  was  greater  than  zero  (stable 
case),  then 

0u  (h/V)    =     1  +  7L^  (67) 

If  L^   was  less  than  or  equal  to  zero  (unstable  case),  then  a  new 
Richardson  number  was  determined  from 

1/4 

Ri  =    L.    (1  -  18  Rj    .  ,)  (68) 

Jnew  *  1old  v 

where  initially  Rj         =   L^  .    This  iteration  was  continued  until  the 

inequality 


xnew        1old 


<      e 


(69) 


was  satisfied,  where     e  =  0.01  for  this  model.    Then  0    (h/V)  was 
found  from 

Vh/L'>   =  71^ TV4  •  (70) 

■^new 
2 .         The  Modified  Blackadar  Baroclinic  Boundary  Layer  Model 
The  basis  of  this  model  was  an  iteration  procedure  which 
included  the  effects  of  stability  and     baroclinicity  on  the  computed 
wind.    The  procedure  was  free  of  computational  instability  and  converged 
rapidly  for  all  reasonable  choices  of  initial  U *  .      First  values   of  zQ  , 
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P/P' /S,q,    a,     7  ,      *V  Q  ,    and  U^  were  obtained  using  respec- 

new 

tively  equations  (23),   (35),  (38),   (43),   (44),   (45),   (46),   (50),  and  (51). 
The  iterative  procedure  was  repeated  until  the  inequality 


U*  -  U* 

new 


<     e 


(71) 


was  satisfied  where    «   was  0.05.    The  iteration  was  repeated,  setting 
,  until  the  inequality  was  again  satisfied.    If  L^  was   not 


u*=u* 


new 


zero,  then  new  values  of     UJ  ,    z_  ,  and  L1 


argument  of     W  was   za/L' .    If 


new 


were  computed  where  the 


L'  -  L' 

new 


>      5 


(72) 


then  L'  was  set  equal  to  L'new  .        The  model  was  repeated  from  where 
L*  was  first  calculated  until  the  inequality  was  satisfied.    After  com- 
pleting the  iterations  for  the  case  where  the  thermal  wind    (r)   and    n 
were  equal  to  zero,  the  computed  values  of  r  and  r|    were  substituted 
and  the  model  was  again  repeated  until  all  criteria  were  satisfied.    If 
IV  0,  then  jef(z/L')  was  computed  as  0(19.5/L')  ;    otherwise  0  (19  .  5/L') 
was  set  equal  to  zero.    Finally  the  wind  at  19.5  meters  was  determined 


from 


U 


u, 


19.5 


In    19,1     .     ^(19^ 
zQ  L1 


(73) 


Values  of  9    -9     ,  G,  z     ,  z/V  ,     i*J     ,    U^  ,  U  log  R    ,  U/G, 

Q.  S  O  O  1  y  »  0  O 

U  /G  and   R|    were  then  available  for  each  grid  point. 
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C.         COMPUTATION  OF  INDEX  OF  REFRACTION  CONSTANT 

In  this  section  possible  relationships  are  considered  between  the 
refractive  index  structure  constant  and  quantities  dependent  only  on  the 
wind  speed  and  air-sea  temperature  difference.    Specifically    those 
quantities  that  can  be  predicted  from  the  model  considered  in  this  study 
are  examined. 

The  propagation  of  light  waves  and  radio  waves  is  affected  pri- 
marily by  air-density  inhomogeneities  and  therefore  by  temperature  fluc- 
tuations.     One  effect  of  these  fluctuations  is  scintillation  or  rapid  changes 
in  intensity  at  a  point.    These  are  distortions  caused  by  variations  in  the 
speed  of  propagation  through  turbulent  elements  having  different  densi- 
ties,  and  hence  different  indexes  of  refraction.    The  refractive  index 
structure  constant  (C    )  is  the  primary  parameter  characterizing  the  turbu- 
lent fluctuations  of  the  atmospheric  refractive  index  for  scales  from  a 

2 
few  millimeters  to  several  meters.    If  the  quantity  Cn    is  known,  cal- 
culations can  be  made  on  scintillation,  beam  spreading,  and  beam  wander, 
all  of  which  are  important  in  describing  light  wave  propagation. 

A  relation  necessary  in  the  formulation  is  the  one-dimensional 
spectrum  for  temperature  fluctuations  in  the  inertial  subrange  described 

by  Corrsin  (1951) 

-1/3  _5/ 

0       (K)  =    B    C  N   K  (75) 

y 

where    c    is  the  rate  of  dissipation  of  turbulent  kinetic  energy,    N  is  a 

measure  of  the  "smearing"  of  temperature  inhomogeneities,    B  is  an 
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empirical  constant  (=0.81,  Boston,   1970),  and  K  is  the  streamwise 

component  of  wave  number.    The  structure  constant  for  temperature  fluc- 

2 

tuations    (CT )  is  defined  as 

2  "V 

CT    =    B     e       '3    N  (76) 

and  is  related  to  the  index  of  refraction  as  follows 

C2    _      jzZtajfil*,*    C2    _  •  (?7) 

a 
Another  expression  which  has  to  be  considered  is  the  nondimen- 
sional  balance  equation  for  turbulent  kinetic  energy,  which,  for  station- 
ary, horizontally  homogeneous  flow  is  (Busch  and  Panofsky,   1968) 

0U    ~    Z/L    "    0t      -    0D   =    0.  (78) 

0     is  the  nondimensional  wind  shear  (8),    0-^  represents  the  diver- 
gence terms ,  and 

*t       =  „3Z  •  <79> 

* 

The  following  relations  arise  from  the  definitions  in  the  non- 
dimensional  turbulent  kinetic  energy  balance  and  definition  of  N 

c  ■  K    -^  (80) 

_  2 
dz  H        a  z 

Substitution  of  these  two  expressions,  plus  the  definition  of  K     (29) 

2 
into  the  definition  of  C      yields  (Panofsky,   1968) 

4  =  »»73  <Jf)^<-'/3<H>2- 
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Equation  (82)  is  not  a  suitable  expression  because  of  the  lack  of 
information  on  0Z       .         However  if  the  effect  of  the  divergence  terms 

is  neglected,  0         can  be  expressed  in  terms  of  0     and  z/L  as 
e  u 

0r      =    0      -    ZA  .  (83) 

c  u 

It  is  noted  that  with  wind-wave  interaction,  the  assumption  that  the 

divergence  terms  can  be  neglected  is  probably  not  valid,  but  the  possible 

2 
effect  this  assumption  has  on  predictions  of  C      has  not  been  studied. 


Substitution    of  (83)  into  (82)  yields 


,-V3 
CT    =    B 


[*u  "  ZA] 


KH        1 


m        u 


4/3    ae  2 

(kz)         ("P").  (84) 
o  z 


This  final  equation  for  C2   lends  itself  to  the  boundary  layer  model  be- 
cause all  the  parameters  can  be  computed  within  the  model.    Thus,  the 

2 
model  can  be  used  to  obtain  a  horizontal  distribution  of  C     over  water 

for  various  wind  and  stability  conditions.    The  average  stability  dQ/dz 
is  determined  from  (14). 

The  ratio  of  the  exchange  coefficients  KTJ/K      is  set  equal  to  1 . 0 

in      m 

in  Cardone's  formulation.    However  it  has  been  shown  that  K_^  and  K 
are  not  usually  identical  except  for  small  values  of    39/  bz    for  which 
the  turbulent  motion  is  not  significantly  influenced  by  buoyancy  forces. 
Businger  et  al.   (1971)  suggested  the  following  form  for  the  ratio  of  eddy 
diffusions 


1/ 
(1-9  ZA)  /2 


a      =     1.35      ^ —T^y       •  (85) 

h  (1-15  z/L)  V4 
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Here  a    =  1.35  for  neutral  conditions,  which  agrees  closely  with  labora- 
tory measurements  (cf.  ,  Hinze,   1959).    A  value  of  1.35  rather  than  1.0 

for  a,    under  neutral  conditions  arises  because  the  von  Karman  constant 
h 

was  equal  to  0.35  instead  of  the  usual  0.40  (Businger  et  al. ,   1971). 


D.        WAVE  HEIGHT  ANALYSIS 

Wave  analysis  is  another  application  for  the  boundary  layer  model. 
For  this  purpose,  a  simplified  version  of  the  FNWC  sea-swell  model 
(Hubert  and  Mendenhall,   1970)  was  utilized.    The  u  and  v  components 
of  the  19 .  5  meter  wind  determined  from  the  boundary  layer  model  were 
used  as  both  current  and  history  input  in  this  analysis. 

The  FNWC  sea-swell  model  computes  a  weighted  wind  from  the 
wind  components  as  follows: 

u    =    0.3  u  +    0.7  ut 

v    =    0.3vt_12    +    0.7vt  (86) 

where  t-12  is  the  12-hour  old  history  wind  speed  and    t  is  the  current 
wind  speed.    These  component  values  are  saved  as  history  for  the  next 
analysis  run,  and  then  the  total  wind    (V)  is  computed.    The  wave  heights 
(W^ )  in  feet  are  then  computed  from 

+    1.0  (87) 


v 


Ui 


Although  this  is  a  very  simple  empirical  method  for  determining  wave 
heights,  the  basic  method  has  been  used  by  FNWC  since  the  mid-1960's 
and  this  particular  model  has  been  used  operationally  since  August  1970. 
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In  addition,  the  FNWC  model  computed  corrections  for  fetch  and 
ice,  and  performed  some  smoothing  of  the  resulting  height  field.    It  was 
assumed  these  corrections  would  be  minimal  for  this  application. 
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IV.      SPECIFICATION  OF  EXTERNAL  PARAMETERS 

A  primary  objective  of  this  study  was  to  apply  the  two-level  baro- 
clinic  boundary  layer  model  to  actual  atmospheric  conditions  and  compare 
the  computed  results  with  the  observed  results.      The  time  period  00  GMT 
28  November  69  through  00  GMT  02  December  69,  with  an  analysis  inter- 
val of  every  12  hours,  was  chosen  for  this  comparison  for  several  reasons. 
First  an  intense  low  pressure  system,  which  influenced  most  of  the  central 
North  Pacific  Ocean,  developed  during  this  period  with  the  system  reach- 
ing maximum  intensity  about  00  GMT  30  November  69.    The  system  pro- 
vided input  data  which  allowed  the  model  to  predict  boundary  layer 
parameters  under  varying  atmospheric  conditions.     Secondly  considerable 
synoptic  data  had  been  gathered  by  the  author  during  an  earlier  study  on 
the  extreme  sea  conditions  (unpublished)  associated  with  this  storm.    The 
latter  gave  a  good  basis  for  comparing  computed  values  to  observed  data. 
Finally,  this  storm  produced  waves  and  swell  which  caused  considerable 
damage  to  several  islands  of  the  Pacific  as  well  as  to  the  spar  buoy  FLIP. 
Therefore  a  wave  generation  model  was  applied  using  the  predicted  sur- 
face winds,  and  these  predicted  waves  were  compared  to  the  observed 
wave  conditions  . 

The   00  GMT   30  November  69  period  was  of  the  most  interest  due 
to  its  extreme  conditions,  extensive  data  coverage,  and  it  was  the  prob- 
able period  when  the  destructive  waves  were  being  generated.    The  various 
corrections  described  above  were  incorporated  into  the  model  for  this 
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time  period,  and  synoptic  data  were  plotted  and  analyzed  for  surface 
pressure,  air-sea  temperature  difference,  surface  wind  speed  and  wave 
height. 

Surface  layer  parameters  were  computed  including  various  combina- 
tions of  stability,  curvature,  and  thermal  wind  to  yield  information  on 
the  relative  importance  of  each  parameter  in  the  formulation.    The  follow- 
ing notation  will  be  used  in  describing  the  results: 

NN    -     non-neutral,  correction  for  stability 

N       -     neutral  -  no  stability  correction 

NC    -     no  curvature  correction  applied 

C       -     curvature  correction  applied 

NT     -     no  thermal  wind  correction  applied 

T        -     thermal  wind  correction  applied 
It  was  recognized  at  the  onset  of  the  study  that  the  specification 
of  the  external  fields  (  AO  and  G)   had  to  be  reasonably  accurate  in 
order  to  evaluate  the  model  on  the  basis  of  observed  wind  fields.     For 
example,  if  the  observed  surface  winds  were  quite  different  than  those 
predicted  by  the  model,  some  estimate  had  to  be  available  on  the  possi- 
bility that  the  discrepancies  were  due  to  poor  specification  of  the  geo- 
stophic  wind  or  the  air-sea  temperature  difference.    Therefore  in  this 
section  consideration  is  given  to  the  relevant  aspects  of  representing 
those  fields  which  make  up  the  external  parameters. 

The  charts  in  the  following  figures  correspond  to  a  1:30  million 
polar  stereographic  projection  with  the  grid  as  depicted  in  Figure  9. 
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A.  FNWC  ANALYSIS 

Sea  level  pressure  and  air-sea  temperature  difference  fields  for 
00  GMT   30  November  appear  in  Figure  10.    A  967  mb  low  pressure 
center  was  located  near  45 °N  180°     with  considerable  thermal  instability 
both  ahead  of  and  behind  the  system.    A  significant  feature  of  this  sys- 
tem was  the  intense  isobaric  gradient  in  the  western  quadrant. 

The  first  consideration  was  the  comparison  of  the  computed  surface 
geostrophic  wind  with  the  computed  gradient  wind  for  this  period.    The 
geostrophic  and  gradient  wind  fields  computed  from  the  surface  pressure 
fields  in  Figure  10  appear  in  Figures   11  and  12  (units  of  meters/second). 
Figure  13  depicts  the  influence  of  the  curvature  correction  on  the  geo- 
strophic wind  computed  as  a  percent  difference  from  geostrophic  wind 
minus  gradient  wind  then  divided  by  geostrophic  wind. 

The  two  maximum  wind  centers  west  and  north  of  the  Low  are  based 
on  data  at  one  grid  point  and  result  from  the  influence  of  anticyclonic 
curvature,   strong  isobaric  gradients,  and  finite  difference  computations. 
These  features  do  not  appear  to  be  realistic.     However,  the  wind  speed 
distribution  depicted  by  them  is  good  except  that  the  magnitude  of  the 
wind  speed  may  be  a  little  too  large.    It  was  expected  that  the  maximum 
change  would  occur  in  an  area  of  maximum  curvature  and  maximum  pres- 
sure gradient  from  the  relation  in   (60).     Generally  speaking,  the  change 
in  wind  speed  due  to  the  curvature  correction  is  10  to  20%  with  a  maxi- 
mum change  of  34%  imposed  by  the  limits  placed  on  the  term  G/Rf  in  (60). 
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B.         SUBJECTIVE  ANALYSIS 

The  results  from  the  hand  analyzed  surface  chart  were  also  used 
in  comparisons  of  results.    Figure  14  illustrates  the  surface  pressure 
and  air-sea  temperature  difference  fields  from  the  subjective  analysis 
which  are  analogous  to  Figure  10.      The  major  variations  between  Figure 
10  and  14  appear  to  be  due  to  the  difference  in  the  number  of  synoptic 
reports  used  as  well  as  the  analysis  scheme.    Theoretically,  Figure  14 
should  be  more  representative  since  more  data  were  included,  but  it  was 
a  subjective  analysis.     Figures   15  and  16  depict  the  geostrophic  and 
gradient  surface  winds  computed  from  the  subjectively  analyzed  pressure 
field . 
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V.      RESULTS  FROM  EXPERIMENTS  WITH  THE  MODEL 

A.  COMPARISONS  IN  TERMS  OF  SURFACE  ROSSBY  NUMBER  AND 

STABILITY  CATEGORIES 

Near  surface  parameters  were  computed  for  the  19.5  meter  level. 
This  level  represents  an  average  height  of  ship  anemometers  (Bunting, 
1968).     Theoretically  the  computed  19.5  meter  level  winds  could   be 
compared  with  the  observed  wind  fields  and  this  would  yield  some 
information  on  the  accuracy  of  the  model  in  predicting  surface  winds. 

In  order  to  verify  the  computations  and  also  to  examine  the  range 

of  conditions  provided  by  the  data,  a  comparison  was  made  with  Cardone's 

results  which  appear  in  Figures  4  and  5.    Results  for  this  comparison 

were  obtained  from  computations  in  which  stability  was  included  but  not 

baroclinicity.    Results  from  all  nine  time  periods  appear  in  Figures  17 

and  18  where  values  of  the  logarithm  of  the  surface  Rossby  number  (24) 

range  from  7.76  to  9.45.    The  U./G  versus  log  R     results  (Figure  17) 

x  o 

agree  with  Cardone's  results  (Figure  4)  except  with  regards  to  the  inter- 
val between  L^  isolines.    In  Figure  17   this  interval  (spacing)  is  larger. 
With  respect  to     ^Q   versus   log  R     (Figure  18),  the  unstable  results 
(negative   L*)    agree  with  those  in  Figure  5  except,  again,  with  regards 
to  the  interval  between   L*  isolines.    The  stable  results  (positive   L^. ) , 
however,  deviate  considerably  from  Cardone's  results  with  the  veering 
angle  decreasing  as    L^  increases.      The  results   in  Figure    18   would 
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indicate  that  the  veering  angle  decreases  with  both  increasing  stability 
and  increasing  instability  with  a  discontinuity  occurring  near  neutral 
conditions . 

The  unusual,  and  unexpected,  features  in  the  results  just  de- 
scribed were  examined  in  considerable  detail.    The  same  results  were 
noted  by  Pafias  (unpublished  study,   1971)  in  an  investigation  with 
Cardone's  program.    These  features  in  the  results  are  perhaps  related 
to  the  drag  coefficient  representation  selected  by  Cardone  and  also  be- 
cause different  iteration  limits  were  used  in  this  study  than  those  used 
by  Cardone. 

The  range  of  data  with  respect  to  U./G  and  9-9     was  also 

*  as 

examined  and  the  results  appear  in  Figure  19.    The  computations  for 
these  results  were  made  on  the  basis  of  a  barotropic  atmosphere  with  no 
correction  for  isobaric  curvature.    These  computations  differ  from  those 
used  by  Cardone  for  Figure  8  in  that  a  variable  coriolis  parameter  was 
used.    Isolines  have  been  drawn  for  geostrophic  wind  speeds  ranging 
from  15  to  75  knots.    For  those  wind  speeds  considered  by  Cardone  in 
Figure  19,  and  also  for  wind  speeds  above  20  knots  considered  by 
Pafias,  these  results  are  in  agreement.    The  change  in  shape  of  the 
curves  for  winds  under  2  0  knots  was  not  described  by  Cardone  but  is 
probably  due  to  a  combination  of  the  functional  forms  for  the  stability 
dependence  and  the  drag  coefficient  as  chosen  by  him. 
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Fig.  19.  G-eostrophic  drag  coefficient  versus  air-sea 

temperature  difference  for  range  of  geostrophic  wind 
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B.         EXPERIMENTS  WITH  CHANGING  THE  PHYSICS  IN  THE  TWO- LAYER 

MODEL  AND  CHANGING  THE  SPECIFICATION  OF  THE  SURFACE 

GEOSTROPHIC  WIND 

In  this  section    U  results  which  were  obtained  from  different 

computations  with  the  two-layer  model  are  compared.  These  computations 
were  varied  by  including  or  not  including  specific  parameters  or  correc- 
tions.    The  variations  considered  and  compared  were  described  earlier 
and  are    1)  N/NC/NT,    2)  NN/NC/NT,    3)    NN/NC/T,    4)  NN/C/T,    and 
5)  NN/NC/T  with  k=0.35. 

These  variations  will  be  referred  to  as  Cases    1  through   5.    It  is 
noted  that  if  the  geostrophic  wind  speed  was  less  than  or  equal  to  5.0 
meters/second,    stability  and  baroclinicity  were  not  taken  into  account 
in  the  iteration  schemes.    In  regions  with  this  wind  speed  range,  espe- 
cially in  the  subtropics  ,  differences  will  not  arise  except  where  isobaric 
curvature  exists . 

Per  cent  difference  (1-2)  will  be  used  to  compare  the  two  cases 
where  (1-2)  corresponds  to  the  following  computation  at  each  grid  point: 
case  1  result  -   case   2  result  nn 

.  A     1UU 

case    1  result 
Although  large  differences  between  cases  were  unexpected,  small  cor- 
rections at  low  wind  speeds  could  produce  large  per  cent  changes. 

The  surface  wind  fields  corresponding  to  each  of  the  five  cases 
are  presented  in  the  following  paragraphs  along  with  selected  comparisons 
between  pairs  of  cases. 
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Case  1  (N/NC/NT) 

This  is  the  simplest  case  with  respect  to  the  possible  physics 
which  are  included  in  the  model.    The  external  wind  field  is  defined  for 
a  neutral,  barotropic  atmosphere  with  no  curvature  correction.    The  wind 
field  at  19.5  meters  for  this  case  appears  in  Figure  20. 

Case   2  (NN/NC/NT) 

The  effects  of  stability  are  now  included  and  the  19.5  meter  wind 

field  for  this  case  appears  in  Figure  21.    The  per  cent  difference  (2-1) 

field  is  shown  in  Figure  22.     Positive  values  correspond  to  increases  in 

the  computed  U        _   field  due  to  including  the  stability  effect.    This  in- 
ly .  o 

crease  is  quite  evident  on  the  western  side  of  the  cyclone.    The  patterns 
are  not  expected  to  be  coincident  with  the  air-water  temperature  isolines 
in  Figure  10  because  the  stability  effect  is  dependent  on  both  the  temper- 
ature gradient  and  the  wind  speed.    The  20%  contours  near  the  eastern 
border  are  attributed  to  small  corrections  at  low  wind  speeds  as  dis- 
cussed previously. 

Case   3  (NN/NC/T) 

The  effects  of  including  the  baroclinicity  or  thermal  wind  are  now 
included  and  the  19.5  meter  wind  field  for  this  case  appears  in  Figure 
23.    The  per  cent  difference  (2-3)  field  is  seen  in  Figure  24.     Negative 
values  correspond  to  an  increase  in  the  computed  19.5  meter  wind  due 
to  including  the  thermal  wind  in  the  model.    In  this  difference  field,  only 
20  grid  points  (6%)  appeared  to  be  significantly  affected  by  including  the 
thermal  wind.    These  grid  points  generally  lie  in  an  area  of  warm  air 
advection.  cc 
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Case  4  (NN/C/T) 

In  addition  to  all  specifications  possible  within  the  physics  of  the 
model,  the  correction  for  curvature  is  applied  to  the  external  wind  field. 
The  19.5  meter  wind  analysis  for  this  field  is  shown  in  Figure  25  and 
the  per  cent  difference  field  (3-4)  appears  in  Figure  26.    Negative  values 
correspond  to  larger  19.5  meter  winds  due  to  including  the  curvature  cor- 
rection.    Strong  wind  maximums,  noted  in  comparing  the  external  wind 
fields  with  and  without  this  correction,  also  appear  in  Figure  26  as  the 
predominant  feature.    In  general,  the  effective  change  in  the  19.5  meter 
wind  field  specification  is  larger  and  covers  more  area  than  any  of  the 
other  cases . 

Case   5  (NN/NC/T  with  k=0.35) 

This  case  differs  from  case  3  only  with  respect  to  the  change  of 
the  von  Karman  constant  from  0.40  to  0.35.    The  19.5  meter  wind  field 
for  this  case  appears  in  Figure  27  and  the  per  cent  difference  field  (3-5) 
is  depicted  in  Figure  28.    Negative  values  correspond  to  increased  19.5 
meter  winds  due  to  the  change  in  the  value  of  k.      The  suggested  value 
of  k=0.35  arose  from  recent  surface  layer  observational  investigations 
by  Businger  et  al.  (1971).    They  noted  from  their  results  that  using 
k=0.40  in  several  empirical  expressions  led  to  estimates  of  U^   values 
which  were  15%  too  high  when  compared  with  direct  measurements  of 
the  momentum  transfer.    The  fields  in  Figure  28  indicate  that  a  general 
increase  of  12  to  17%  occurs  in  the  computed  19.5  meter  wind  field  if 
the  lower  value  of  k  is   used.    This  increase  agrees    with   the  U^ 
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discrepancies  noted  by  Businger  et  al. .      Because  the  von  Karman 
constant  appears  in  several  expressions  in  the  boundary  layer  model, 
the  importance  of  its  change  in  the  two-layer  model  was  not  known. 
These  results  indicate  that  the  difference  between  k=  0.35  and  k=  0.40 
is  important. 

In  summary,  the  five  cases  described  above  indicate  which  terms 
were  most  influential  in  the  model.    It  should  be  remembered  that  most 
of  the  subtropical  points  were  not  influenced  by  the  model  since  the 
wind  speeds  were  generally  less  than  5.0  meters/second.    Of  all  the 
modifications  to  the  model,  the  curvature  correction  to  the  external 
field  produced  the  largest  variation  over  the  entire  grid.     Changing   k 
from  0.40  to  0.35  accounted  for  the  second  most  significant  change. 
Differences  in  the  surface  wind  fields  were  evident,  however,  when 
specifications  in  the  model  were  made.     Of  the  two  possible  specifica- 
tions ,  stability  appeared  to  be  more  significant  than  the  thermal  wind 
with  respect  to  area  influenced  and  per  cent  difference  changes . 

Although  the  influence  of  the  thermal  wind  in  this  model  appears 
to  be  small  with  respect  to  the  relative  magnitude  of  the  computed  wind, 
the  influence  it  has  on  the  direction  of  the  wind  appears  to  be  signifi- 
cant.   The  difference  field,  in  degrees,  of  the  wind  direction  for 
NN/NC/T  -  NN/NC/NT  appears  in  Figure  29.    This  field  represents  the 
difference  in  wind  direction  between  a  case  with  and  a  case  without  the 
thermal  wind  specification.    The  region  of  significant  differences  does 
not  coincide  with  that  region  where  the  thermal  wind  effect  was    large 
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with  respect  to  wind  speed.    Positive  values  in  Figure  29  correspond  to 
a  veering  wind  and  negative  values  correspond  to  a  backing  wind.    It  is 
noted  that  the  direction  change  due  to  baroclinicity  is  primarily  deter- 
mined by  the  angle  between  the  thermal  and  geostrophic  winds. 

C.         COMPARISONS  OF  COMPUTED  WINDS  WITH  OBSERVED  WINDS 

These  comparisons  will  be  described  for  the  case  NN/C/T  with 
k=0.35.    The  19.5  meter  wind  fields  computed  from  the  other  cases  were 
compared  to  the  subjectively  analyzed  observed  surface  wind  field 
(Figure  30)  and  in  all  cases  the  computed  winds  underestimate  the  ob- 
served winds.    This  was  also  the  case  when  the  above  corrections  were 
included,  but  yielded  the  best  results  with  respect  to  the  magnitude  and 
patterns  of  the  wind  speeds. 

Because  the  previous  comparisons  revealed  the  importance  of  the 
specifications  of  the  external  wind  field  with  respect  to  curvature,  it 
was  considered  necessary  in  this  comparison  to  examine  the  importance 
of  the  initial  analysis  (FNWC  or  subjective).      The  19.5  meter  wind  field 
computed  with  the  above  specifications  and  with  the  FNWC  analysis  is 
shown  in  Figure  31  and  the  wind  field  based  on  subjective  analysis  is 
seen  in  Figure  32.    The  observed  surface  wind  field  which  appears  in 
Figure  30  was  analyzed  using  only  those  observations  reported  above 
10  meters/second.    In  Figures  31  and  32,  the  maximum  contour,  30 
meters/second,  agrees  with  the  observed  maximum  wind  speed,  but  the 
objective  analysis  does  not  indicate  the  25  meter/second  maximum  east 

of  the  low . 
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A  comparison  of  computed  19.5  meter  winds  from  the  two  analyses 
(Figures  31  and  32)  and  the  external  wind  fields    (Figures  11  and  15) 
shows  that  differences  in  the  external  wind  fields  are  consistently 
carried  through  the  model.    The  model,  however,  was  able  to  produce 
19.5  meter  wind  fields  in  both  cases  which  were  in  good  agreement  with 
the  observed  wind  field,  at  least  for  wind  speeds  over  10  meters/second 

D.         RESULTS  FROM  WAVE  HEIGHT  COMPUTATIONS 

Wave  height  fields  were  computed  using  the  19.5  meter  wind 
fields  obtained  from  the  boundary  layer  model  and  the  simplified  wave 
height  program  developed  from  the  FNWC  operational  sea-swell  model. 
The  program  was  run  starting  at  00  GMT  28  Nov  in  order  to  build  history 
for  this  analysis.    The  19.5  meter  wind  computations  were  obtained 
using  all  specifications  and  corrections  described  previously,  and  they 
were  also  made  on  both  the  FNWC  and  subjectively  analyzed  fields.  The 
computed  wave  field  from  the  FNWC  data  is  seen  in  Figure  33  and  the 
wave  field  from  the  subjective  analysis  appears  in  Figure  34.    The  ob- 
served wave  field  obtained  from  synoptic  observations  is  depicted  in 
Figure  35.    It  is  noted  that  wave  height  is  probably  the  most  difficult 
parameter  to  observe  and  therefore  a  representative  analysis  of  the  ob- 
served wave  field  is  difficult  to  achieve. 

The  wave  field  obtained  from  the  subjective  analysis,  Figure  34, 
appears  to  be  in  good  agreement  with  the  observed  heights  with  maxi- 
mum computed  heights  slightly  less  than  the  maximum  observed  heights. 
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Results  computed  from  the  FNWC  analysis,  Figure  33,  appear  to  have 
similar  agreement  with  respect  to  the  maximum  height  but  poorer  agree- 
ment with  respect  to  the  overall  pattern  of  the  wave  field. 

In  general  these  results  indicate  that  the  two- layer  model  can 
produce  surface  wind  fields  which  can  be  used  for  wave  height  analysis, 
even  from  a  very  simple  wave  model.    It  is  noted  that  in  the  subtropics 
the  wave  heights  computed  from  both  sets  of  analysis  are  considerably 
below  the  observed  wave  heights.    This  is  probably  due  to  low  surface 
geostrophic  wind  speeds  in  this  area  and  the  fact  that  no  corrections  are 
made  for  stability  or  baroclinicity  at  low  wind  speeds.    This  could  pos- 
sibly be  corrected  in  the  model  by  selecting  a  higher  initial  value  for 
the  drag  coefficient. 

E.  RESULTS  FROM  INDEX  OF  REFRACTION  COMPUTATIONS 

A  final  application  of  the  boundary  layer  model  was  the  computa- 
tion of  the  index  of  refraction  structure  constant  (C_?)  using  (84).  The 
ratio  of  the  exchange  coefficients  was  determined  from  (85)  and  the 
input  data  for  this  equation  were  obtained  by  including  the  stability 
effect,  the  thermal  wind  effect,  the  correction  for  curvature,  and  using 
k=0.35.  All  parameters  needed  to  compute  C^  were  determined  in  the 
model. 

The  results  from  these  computations  appear  in  Figure  36.    Values 

2 
of  C     range  from  0 . 0  in  regions  of  near  neutral  stratification  to  a  maxi- 
mum of  0. 181  in  the  region  of  maximum  wind  speeds.    Comparisons   of 
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these  values  with  those  measured  by  Portman  et  al.  (1968)  indicate  that 

they  are  comparable  in  regions  of  low  wind  speed  with   0-0      less  than 

a        s 

2 
1   degree.    Although  no  definite  relation  between  C     and  wind  speed  or 

2 

Q    -Q     appears  in  Figure  36,   larger  values  of  Cm  appear  to  be  associ- 
a       s  T 

ated  with  the  greater  wind  speeds  and  larger  instability  (€    -  9     ranging 

a        s 

from  -1  to  -6°C). 

2 

The  results  indicate  that  a  horizontal  distribution  of  C      for  a 

given  height,  in  this  case    19.5   meters,    can  be  computed   with  the  model 
with  some  degree  of  confidence  in  the  range  of  values  obtained.    Com- 
parisons should  be  made  between  results  computed  by  this  model  and 

2 
observations  of  CT  to  obtain  better  verification. 
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VI.      CONCLUSIONS 

The  results  of  this  study  have  indicated  that  a  marine  boundary 
layer  model  can  be  applied  to  synoptic  scale  data  to  obtain  various  sur- 
face parameters.    These  parameters  include  friction  velocity   (U^.), 
surface  roughness    (z    ),    surface  Rossby  number   (R    ),    19 . 5  meter  wind 
(U  ),    veering  angle    (^„),  and  stability  length   (z/L) . 

First,  the  results  indicate  that  the  curvature  correction  produced 
the  greatest  change  to  the  wind  analysis.    However,   stability  and  modi- 
fication of  the  von  Karman  constant  produced  significant  changes  in  the 
field,  while  changes  due  to  the  thermal  wind  were  confined  primarily  to 
changes  in  wind  direction. 

Second,  the  19.5  meter  wind  was  computed  including  various 
physical  specifications  and  corrections  and  compared  with  synoptic 
observations.    The  results  indicate  that  the  computed  wind  was  lower 
than  the  observed  wind  in  all  areas  and  for  all  cases,  but  the  general 
patterns  agreed  well  with  observations . 

Finally,  it  appears  that  Cardone's  planetary  boundary  model  pro- 
vides a  reasonable  framework  to  specify  those  surface  layer  parameters 
required  for  spectral  wave  forecasting  and  index  of  refraction  calculations. 

Further  studies  are  needed  to  more  completely  verify  and  improve 
the  model.    The  specification  of  the  drag  coefficient  is  probably  the 
aspect  in  this  model  which  needs  the  most  attention.    The  drag  coefficient 
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has  received  considerable  examination  in  observational  investigations 
and  the  results  indicate  some  disagreement  among  investigators.    A 
suggestion  of  this  study  is  to  increase  the  initial  value  of  C     to  increase 
the  computed  surface  wind  speeds  especially  in  areas  of  low  wind  speed. 
Application  of  this  model  to  a  region  of  reliable  synoptic  data  such  as 
during  the   BOMEX   experiment  in  1969   would  reduce  some  of  the  errors 
due  to  data  input  and  also  test  the  physics  of  the  model  in  the  tropical 
atmosphere.     Finally,  observations  of  the  index  of  refraction  over  a  wide 
range  of  wind  speed  and  stability  would  allow  reliable  comparisons  to  be 
made  of  the  structure  constant   (C-p  ) . 
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